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N-body simulations of dark matter halos show that the density profiles of halos behave as p(r) oc r~ a ( r \ 
where the density logarithmic slope a ~ 1 ~ 1.5 in the center and a ~ 3 ~ 4 in the outer parts of halos. 
However, some observations are not in agreement with simulations in the very central region of halos. The 
simulations also show that velocity dispersion anisotropy parameter f3 in the inner part of the halo and 
the so called "pseudo phase-space density" p/er 3 behaves as a power-law in radius r. With these results in 
mind, we study the distribution function and the pseudo phase-space density p/cr 3 of the center of dark 
matter halos and find that they are closely-related. 
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1. Introduction 

The formation and evolution of the dark matter halo which can be treated as the self-gravitational 
collisionless stellar system have become a challenging issue in the study of dark matter. 

Thanks to the improved computational power, N-body simulations of dark matter halos become 
more and more accur ate and important with increasing resolution. N-body simulations such as the 
universal NFW profile^^ and other JSEEI show that the density profiles of dark matter halos behave 
as r~ a ( r \ where a £± 1 ~ 1.5 in the center and a ~ 3 ~ 4 in the outer parts of halos. However ,the 
numerical inner behaviors of dark matter halos are not supported by observations! I I I I I I I I 
Some work indicates that density profiles of dark matter halos might b ecome shallower than r ~ 1 
in the center^, and perhaps even tend to be a core with no cusp at all! ^^ 

N-body simulations not only provide us the density profiles of halos, but also give the rele- 
vant information of the velocity space of collisionless particles in the halos. Velocity dispersion and 
anisotropy profileA ^^^^^^^ nave been well described by simple analytical fits. Two inter- 
esting phenomena from simulations indicate that vel ocity d ispersion and the density profiles of the 
haloes are not independent. First, Hansen and Moor jMEU f oun d that the density logarithmic slope 
a{r) is correlated to the velocity anisotropy which is parameterized by the anisotropy parameter 
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(3(r) and they provided the empirical formula f3 « 1 — 1.15(1 — a/6). Therefore f3 « (isotropic 
velocity dispersion) in the inner part as a ~ 1 and /3 w 0.5 in the outer part as a f=a 3. Second, it has 
been argued^ that the so called pseudo phase-space density follows a power law p(r)/a 3 (r) oc r 1 
with exponent 7 = 1.875, where p(r) is the density profile and a 2 is the total velocity dispersion. 
Subsequent studies have confirmed that p(r)/a 3 ( r) is a powe r law in radius, but the best fitting 
values of the exponent 7 diverse from each 

othei |17 | 27 | 28 | 29 | 

and range from 7 = 1.90 ± 0.05 to 
2.19 ± 0.03. Because p/a 3 has the same dimension as the phase-space density, p/a 3 has been called 

loni 

pseudo phase-space density or "poor-man's" phase-space density^. 

With these results in mind, much theoretical work has been done for the study of the relation 

between density profile behavior and pseudo phase-space density. Some aut hors examine this matter 

by solving the Jeans equation. Williams et al got a critical exponents^ 7 = 35/18 and Dehnen 
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and McLaughlin calculated corresponding density profiles in both isotropic and anisotropic casea^ 1 . 
Some other work solved t he Jeans equ ation to explore the relation between density profile and 
pseudo phase-space density | 14 | 30 | 33 |5H. Recently, R. N. Henrikserf23 considered a series expansion 
for a dark matter distribution function in the spherically symmetric anisotropic limit to discuss 
pseudo phase-space density. 

In this paper, we concentrate on the center of the dark matter halo where the velocity dispersion is 
almost isotropic and calculate the distribution function and pseudo phase-space density where r — > 
in the spherically symmetric case. In Section 2 we review the basic knowledge about distribution 
function which is needed for this paper. In Section 3, the distribution function, velocity dispersion 
and pseudo phase-space density in the center of the dark matter halo are calculated with the given 
density profile. We make the discussion and conclusion in Section 4. 



2. General Formulae 

In the center of dark matter halos, with spherically symmetric assumption, the distribution function 
F(E,L) can be reduced to F(E) as anisotropy parameter /3 w (almost isotropic) when r — ► 0. If 
we know the distribution function, we can calculate the density profile and total velocity dispersion 
profile as below^", 



p(iP) = An / F(E)^/2(^-E)dE , (1) 



a 2 = al + 4 = tllnMu* f F(E)(ij; - Ef/ 2 dE, (2) 
p Jo 

where the binding energy E is defined as 

E = ^{r)-\v 2 r -\v 2 T , (3) 
and tp( r ) is t ne relative gravitational potential which can be obtained from the Poisson equation: 



1 d ( 2 dip 



r ' _L = ^Gp{r) . (4) 
dr \ ar 



v r is the radial velocity and vt is the tangential velocity: 



VT = J v 2 g +vl. (5) 
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EddingtorP3 provided the inversion formula of Eq. ([I]) : 

where D% denotes the 2nd order differentiation operator with respect to E. 
The anisotropy parameteJ^ mentioned in the Introduction is defined as: 

where a\ and of are the tangential and radial velocity dispersion. If f3 < 0, the velocity dispersion 
is tangentially anisotropic. < (3 < I correspond to the radially anisotropic cases, and /? = 0, the 
isotropic case. 



3. Distribution Function &: Pseudo Phase-Space Density 

In this Section, we will calculate the distribution function and pseudo phase-space density with 
the given density profile in the very central region of the spherically symmetric dark matter halo 
where the velocity dispersion is almost isotropic. Let us consider a family of density profiles with 
parameters a and e named as New generalized NFW profiles 

P ~ C (r/r s ) a (l + r/r s ) 3+£ -° ^ 

We set the characteristic radius r s — 1, the total mass Mtot = 1 and the gravitational constant 
G = 1 here. Then the density profiles reduce to 

P = C—, ^-rr , (9) 

Then the relative potential ip(r) can be calculated from the Poisson Eq. 

V(r) = ^0 - 4nCr- a r(-a) (A^r) - 2A 2 (r) + A 3 (r)) 

A 1 (r)= 2 Fi[l-a + e, -a, 2-a, -r] (11) 
A 2 {r) = 2 Fi[2-a + e, -a, 2-a, -r] 
A 3 {r)= 2 Fi[3-a + e, -a, 2-a, -r] 

where ipo = 2Z~~ ^ s the relative gravitational potential at r = which is determined by the condition 
V , ( r )l r ^ o = 0- ^nd then, we calculate the asymptotic approximation of VKOIr- >o anc ^ P(V0lr->o as 
below: 



4 D. Ma and P. He 



6-5« + a» ^ (13) 
47rC 

From ip{ r )\ r ^o an d /°('0)lr^O' we can S e t the asymptotic approximation of F(E)\ r _ tQ when using 
Eq. © 

i?< /(a ' 2) [(a - 2)V> (( tt - 2)^o + 4g + (4 - 4a - Za^^y/^) ff 2 A 4 (E)} 
( )lr-.o- 4\/2(2-a) 2 £ 3 / 2 (£;-^o) 2 7r 2 



, , x r l a 3 E 

ME) = 2F 1 [-, 



2' 2- a 2' Vo 

(14) 

We should note that this result of distribution function is only valid in the very center of dark 
matter halo where the velocity dispersion is almost isotropic. In the outer part of the halo, no 
matter whether the velocity dispersion is isotropic or not, this result is invalid and even unphysical. 
It's easy to note that, when r — > 0, the binding energy E approaches to its maximum value i/'o- 
Following this observation, we can further reduce Eq. (|14[) in the parameter space | < a < 2 which 
has covered the result of simulations. Nevertheless, it's just in this parameter scope that F(E) in 
Eq. ([T4"|) is physically meaningful: 

2 
3 

By using equations ([2]), (fi"3]) and (fl5|) . we can get the total velocity dispersion in the very center 
of the halo, 



F(E)\ e ^ q = d(V>o - E)-*-^, - < a < 2 (15) 



2/ / \ I 2 9 / 2 dW 5/2 M M ^_J>- Fr 5 1 2 7 V 



^-^o 5B vru r; L 2' 2 o-2' 2' V-V'o 



r5 1 2_ 7 j> ■ 

1.2' 2 a-2' 2' i/>o- 



_ 5(a-2) , g \-H^5 ( 16 ) 
^0 4(a-l)^ -V 

15y / 7rr(-2-^) , g \-5/2 
1 a_1 U/, n _,/ J 



Then, it's not difficult to simplify equation (| 16[) to the final result: 

(ifo-il)), 2>o>l, 



<^(voi^ oc < (l?) 



(V>o - V') 2 ""* 1 > a > 2 



3 ' 



and substitute equation (| 1 2|) into equation (|17jl to transform <r( , 0) into cr(r) 



r 2 " Q , 2>a>l, 
^WUn^ (18) 



r a , l>a>|. 
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Now, after we know the total velocity dispersion and the density profile in the very center of the 
halo, the pseudo phase-space density can be calculated directly, 



p{r) 



—f 




oc r 




2 > a > 1, 



(19) 
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1 > a > §. 



The power-law behaviour of the ps eudo phase- space density was first found by simulation^ and 
analyzed by subsequent theoretic wor l|201 32 | 33 | 34 |^ its physical meaning is still unclear. For this 
consideration, let us compare the distribution function and pseudo phase-space density in the very 
center of the halo. 

We should note that the velocity of particles, which may stably exist in the very center of the 
halo, approaches to with r — > 0. If the velocity term of the binding energy Eq. or the total 
velocity dispersion a 2 is much smaller than the second term of the RHS of Eq. (TT2"]) or at most has 
the same asymptotic behavior as the second term of RHS of Eq. ([12")) , this condition requires a > 1 
which can be obtained from Eq. (fT2"f and Eq. (fTS)) , then (ipQ — E) \ r ~,o oc r 2 ~ a . With this formula 
and Eq. (fl~5*|) , we can get 



From Eq. (fTO)) and Eq. (f2"D|) , we see that, if 2 > a > 1, the distribution function and the pseudo 
phase-space density have the same asymptotic behavior where r — > 0. This result indicates that the 
"poor-man's" phase-space density is not so poor but might relate to real phase-space density at least 
in the very center of the halo. 

4. Discussion and Conclusions 

N-body simulations show us some properties of the dark matter halo, such as density profile p(r), 
anisotropy parameter (3(r), pseudo phase-space density, and so on. Although the simulation's reso- 
lution is limited, these information should be considered properly in the theoretic analysis in order 
to make theoretic work more realistic. 

In this paper, we investigate the problem of the cente r of the halo based on the following facts: 
(1) the existence of universal density profiles such as NfWEEI Moored, and in particular the New 
generalized NFW profile used in our work; (2) anisotropy parameter (3 s» in the center of the halo; 
(3) the pseudo phase space density follows a power law in radius r; (4) limited resolution in the 
central region. 

We should notice that in this paper all results are only valid in the very center of the halo. From 
a given density profile, we calculate the asymptotic approximation of the distribution function in 
the very center of the halo. Then the total velocity dispersion and the pseudo phase-space density 
are obtained from the distribution function. Eq. (fl9|) shows the relation between the two parameters 
7 and a. If we set 7 = 1.94 as some simulations indicate, then a should be 0.776 which is smaller 
than NFW's result. Otherwise, if we set a ~ 1 ~ 1.5 to adapt NFW profile and Moore's profile, then 
7 ~ 2.25 ~ 2.5 which is larger than the simulations' result. This apparent contradiction between 



F(E)\ 



oc r 



-3+f 



2, a > 1. 



(20) 
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theoretic result and simulation is not so weird since the resolution in the central region isn't high 
enough to describe the very central region of the dark matter halo. 

Comparing the distribution function (real phase-space density) with the pseudo phase-space 
density, we find that they have the same asymptotic behavior if 2 > a > 1. This result indicates 
that the distribution function and pseudo phase-space density are closely related, even though the 
assertion is rather premature that the distribution function might have a power-law behavior just 
as the pseudo phase-space density in the whole range of radius r. This interesting suggestion may 
also give us more confidence in the study of pseudo phase-space density and some new clues to 
the construction of the distribution function. By the way, some authors^ a l so study the pseudo 
phase-space density in from p/a 3n . It is obvious that our conclusion in this paper is well-founded 
only if n = 1. So, from this point of view, p/er 3 may be the best choice of the pseudo phase-space 
density. 

For the intensive study of the dark matter halo, the full range of radius r should be involved. We 
need to consider how to extend the relationship between "real" and "pseudo" phase-space density 
from the very central region to other regions where the velocity dispersion is anisotropic. Further- 
more, reasonable knowledge and new strategies should be introduced and developed for the con- 
struction of more realistic models. First, it is admitted that dark matter halos in the real world arc 
not spherically symmetric. Second, the physics in the center of the halo may be very complicated. 
Third, the halo actually is a polycomponent system which contains the dark matter, baryonic matter 
and even a supermassive black hole in the center. It is not radical to say that new methods and 
ideas are still needed in the future research on this matter. 
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